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Abstract. I discuss recent calculations of the thermal suppression of bottomonium states in rel- 
ativistic heavy ion collisions. I present results for the inclusive Y(ls) and T(2s) suppression as a 
function of centrality. I compare with recent CMS preliminary data available at central rapidities 
and make predictions at forward rapidities which are within the acceptance of the ALICE dimuon 
spectrometer. 
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INTRODUCTION 

The purpose of ongoing and upcoming heavy ion collision experiments at the Relativistic 
Heavy Ion Collider (RHIC) and the Large Hadron Collider (LHC) is to study the 
behavior of nuclear matter at high energy density, £>1 GeV/fm . At such high energy 
densities one expects to create a deconfined quark gluon plasma (QGP). One of the 
key observables measured in the experiments is the relative yield of heavy quarkonium 
production relative to that in proton-proton collisions scaled by the number of nucleon 
collisions, Raa - Early studies predicted that at high temperatures Debye screening would 
lead to the dissolution of hadronic bound states [1]. Therefore, one might be able to use 
the relative suppression of heavy quark bound states as a "smoking gun" for the creation 
of the quark gluon plasma. Over the last 25 years most of the interest has been focused 
on the suppression of heavy quark-antiquark bound states following early predictions of 
thermal J suppression [2, 3]. In recent years interest has shifted to bound states of 
bottom and anti-bottom quarks (bottomonium) for the following reasons 

• Bottom quarks (m^ ~ 4.2 GeV) are more massive than charm quarks (m c ~ 
1.3 GeV) and as a result the heavy quark effective theories underpinning phe- 
nomenological applications are on much surer footing. 

• The masses of bottomonium states (my ~ 10 GeV) are much higher than the 
temperatures (T < 1 GeV) generated in relativistic heavy ion collisions. As a result, 
bottomonium production will be dominated by initial hard scatterings. 

• Since bottom quarks and anti-quarks are relatively rare within the plasma, the 
probability for regeneration of bottomonium states through recombination is much 
smaller than for charm quarks. 

• Due to their higher mass, the effects of initial state nuclear suppression are expected 
to be reduced. 

As a result one expects the bottomonium system to be a cleaner probe of the quark gluon 
plasma than the charmonium system for which the modeling has necessarily become 



quite involved. For this reason I will focus on the bottomonium states in this paper and 
only consider the thermal suppression of these states, ignoring initial state effects and 
any possible thermal generation or recombination. 

In this paper I will review recent theoretical calculations of bottomonium suppression 
at energies probed in relativistic heavy ion collisions at the Large Hadron Collider 
(LHC). I will present a brief overview of the important aspects of the calculation and 
refer the reader to Refs. [4, 5] for details. I will compare my prediction for inclusive 
T(ls) and Y(2s) suppression with recent preliminary data to be released from the CMS 
collaboration. In addition, I present new predictions for T(ls) and T(2s) suppression 
at forward rapidities in order facilitate comparison with results forthcoming from the 
ALICE collaboration. 



In recent years there have been important theoretical developments in heavy quarko- 
nium theory. Chief among these are the first-principles calculations of imaginary-valued 
contributions to the heavy quark potential. The first calculation of the leading-order per- 
turbative imaginary part of the potential due to gluonic Landau damping was performed 
by Laine et al. [6]. Since then an additional imaginary- valued contribution to the poten- 
tial coming from singlet to octet transitions has also been computed using the effective 
field theory approach [7]. These imaginary- valued contributions to the potential are re- 
lated to quarkonium decay processes in the plasma. The consequences of such imaginary 
parts on heavy quarkonium spectral functions [8, 9], perturbative thermal widths [6, 10], 
in a T-matrix approach [1 1-13], and in stochastic real-time dynamics [14] have recently 
been studied. 

In addition, there have been significant advances in the dynamical models used to 
simulate plasma evolution. In particular, there has been a concerted effort to understand 
the effects of plasma momentum- space anisotropics generated by the rapid longitudinal 
expansion of the matter along the beamline direction. Dynamical models are now able to 
describe the anisotropic hydrodynamical evolution in full (3+l)-dimensional simulations 
[15-20]. This is important because momentum- space anisotropics can have a significant 
impact on quarkonium suppression since in regions of high momentum- space anisotropy 
one expects reduced quarkonium binding [21-25]. In Refs. [4, 5] the dynamical evolu- 
tion of the anisotropic plasma was combined with the real and imaginary parts of the 
binding energy obtained using modern complex- valued potentials. 



The heavy quark potential has real and imaginary parts, V = 9?[V] + z'3 [V]. Here I will 
focus on a model in which the real part of the potential is obtained from internal energy 
of the system since models based on the free energy seem to be incapable of reproducing 
either the LHC or RHIC Raa \?] ■ The real part of the potential is given by [5] 



THEORETICAL METHODS 



Potential model and binding energies 
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where a = 0.385 and o = 0.223 GeV 2 [26] and the last term is a temperature- and 
spin-independent finite quark mass correction taken from Ref. [27]. In this expression 
fl = 5f , 9)m,D is an anisotropic Debye mass where £f is a rather complicated function 
which depends on the degree of plasma momentum- space anisotropy, t, , and the angle 
of the line connecting the quark-antiquark pair with respect to the beamline direction, 6, 
and m D is the isotropic Debye mass [5]. In the limit q~ — > one has = 1. 

The imaginary part of the potential 3 [V] is obtained from a leading order perturbative 
calculation which was performed in the small anisotropy limit 

3[V] = -a s C F T {(j>(r/m D ) - !g [Wl(r/m D , 0) + y 2 (r/m D , 0)]} , (2) 

where <j), y/i, and ty/2 can be expressed in terms of hypergeometric functions [24]. 

After combining the real and imaginary parts of the potential, the 3d Schrodinger 
equation is solved numerically to obtain the real and imaginary parts of the binding 
energy as a function of t, and phard [25, 28]. The imaginary part of the binding energy is 
related to the width of the state 

n „ f 2Z[E hind (T,x ± ,g)} 9t[£ b ind(T,x ± ,g)]>0 fV . 

HT,x ± ,g- j 1QGeV 9i[£ b i„d(T,x ± ,g)]<0 (i) 

where g = arctanh(z/?) is the spatial rapidity. The value of 10 GeV in the second case is 
chosen to be large in order to quickly suppress states which are fully unbound (which is 
the case when the real part of the binding energy is negative). 



Dynamical Model 

The dynamical model used gives the spatio-temporal evolution of the typical trans- 
verse momentum of the plasma partons, Phard^x), and the plasma momentum- space 
anisotropy, £(t,x), both of which are specified in the local rest frame of the plasma. 
The widths obtained from solution of the 3d Schrodinger equation are then integrated 
and exponentiated to compute the relative number of states remaining at a given proper 
time. This quantity is then averaged over the transverse plane taking into account the 
local conditions in the plasma and weighting by the spatial probability distribution for 
bottomonium production which is given by the number of binary collisions computed in 
the Glauber model with a Woods-Saxon distribution for each nucleus. For the temporal 
integration the initial time is set by the formation time of the state in question. 

The resulting Raa is a function of the transverse momentum, pj, the rapidity g, and the 
nuclear impact parameter b. To compare to experimental results transverse momentum 
cuts are applied using a n(pj) = hqE t 4 spectrum. In addition, any cuts on the rapidity 
due to detector acceptance and centrality are applied as needed. For details of the 
dynamical model and Raa computation I refer the reader to Ref. [5] . 



Initial conditions 



For the initial conditions I use a Woods-Saxon distribution for each nucleus and 
determine the transverse dependence of the initial temperature via the third root of 
the number of participants (wounded nucleons). In the spatial rapidity direction I will 
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FIGURE 1. Predictions for the central rapidity inclusive (a) Y(ls) and (b) Y(2s) suppression including 
feed-down as a function of A^ part along with recent preliminary data from the CMS collaboration. Three 
curves show the variation with the assumed shear viscosity to entropy ratio. 



investigate two possible temperature profiles: (a) a broad plateau containing a boost- 
invariant central region with Gaussian limited-fragmentation at large rapidity 

n(g) = n exp(-(|g| - g flat /2) 2 /2o 2 q ) 0(|g| - g flat /2) , (4) 

where £fl at = 10 is the width of the central rapidity plateau, Og = 0.5 is the width of the 
limited fragmentation tails, and uq is the number density at central rapidity [29]; and (b) 
a Gaussian motivated by low-energy fits to pion spectra 

n(g) =n exp(-g 2 / 2 C7?) with cr 2 = 0.64 • 8c 2 ln(v^/2m p ) /3(1 - 4) , (5) 

where c s = 1/ a/3 is the sound velocity, m p = 0.938 GeV is the proton mass, and y/s NN 
is the nucleon-nucleon center-of-mass energy [5, 30]. The temperature distribution is 
given by T ~ n 1 / 3 . 1 note that (4) has the advantage that it has been tuned to successfully 
describe the rapidity dependence of the elliptic flow in LHC heavy ion collisions. 



RESULTS AND CONCLUSIONS 

When considering the suppression of the T(ls) and T(2s) states it is important to in- 
clude the effect of feed-down from higher excited states. In pp collisions only ap- 
proximately 51% of T(ls) states come from direct production and similarly for the 
T(2s). One can compute the inclusive suppression of a state using ^^[Y^s)] = 
Lie states fiRi,AA where fi are the feed-down fractions and is the direct suppression 
of each state which decays into the T(ns) state being considered. Here I will use /; = 
{0.510,0.107,0.008,0.27,0.105} for the T(ls), J(2s), J(3s), % bl , and %b2 feed-down 
to T(l^), respectively [31]. For the inclusive T(2*) production I use = {0.500,0.500} 
for the T(2*) and T(35) states, respectively. For details of the computation of the direct 
Raa for each state see Ref. [5]. 




FIGURE 2. Predictions for the inclusive (a) Y(ls) and (b) Y(2s) including feed-down as a function of 
Afpart at forward rapidities 2.5 < y < 4.0. Bands show variation with rj /S with the bottom, middle, and top 
lines in each set corresponding to 4nrj/S — {1,2,3}, respectively. 



In Fig. 1 I compare to preliminary data on the inclusive Y(\s) and Y(2s) suppression 
available from the CMS collaboration [32]. For this figure I used a broad rapidity plateau 
as the initial density profile as specified in Eq. (4). The central temperatures were taken 
to be T = {520, 504, 494} MeV at t = 0.3 fm/c for 4^T] /S = {1,2,3}, respectively, in 
order to fix the final charged multiplicity to dN^/dy ~ 1400 in each case. As can been 
seen from this figure, the predictions agree reasonably well with the available data. The 
data seem to prefer the largest value of r\ /S shown; however, there is a ±14%, ±21% 
Is, 2s global uncertainty reported by CMS, making it hard to draw firm conclusions. 

Next I consider forward rapidities. The ALICE forward muon spectrometer detects 
charmonium and bottomonium states via dimuon decay with an acceptance that covers 
the pseudorapidity interval 2.5 < y < 4.0 with a resolution of approximately 100 MeV. 
For more information on ALICE'S dimuon capabilities see e.g. [33-36] and references 
therein. One of the key questions that such forward rapidity measurements will be able 
to answer is the nature of the temperature distribution of the plasma at forward rapidities. 

In order to illustrate this in Fig. 2 I show predictions for the inclusive suppression 
of T(l.s') and T(2s) states using either Eq. (4) or Eq. (5) as the initial rapidity profile 
of the density. As can be seen from this figure there is a rather large effect. The boost- 
invariant plateau with limited fragmentation (4) is the one which is canonically used in 
heavy-ion models and follows naturally from the Bjorken picture, while the Gaussian 
assumption (5) is more in line with the Landau picture. I would expect, based on ability 
of the Bjorken-like picture to describe the rapidity dependence of elliptic flow [29], that 
this is, in fact, what one would see in the experiment; however, for completeness sake I 
show the Gaussian in order to demonstrate that the future ALICE results in this region 
of phase space could easily discern between the two possibilities. 

To conclude, in this paper I have attempted to review recent calculations of bottomo- 
nium suppression in LHC heavy ion collisions. I presented comparisons of prior predic- 
tions with CMS preliminary data for inclusive Y(l,s) and T(2s) suppression. In addition, 
I have made new predictions for inclusive T(ls) and T(2s) at forward rapidities within 



the ALICE acceptance. I found that one can easily distinguish between the so-called 
Landau and Bjorken pictures of heavy ion collisions with such forward rapidity data. 
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